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In this paper, macroscopic quantum tunneling (MQT) effect of Z 2 topological 
order in the Wen-Plaquette model is studied. This kind of MQT is characterized 
by quantum tunneling processes of different virtual quasi-particles moving around a 
torus. By a high-order degenerate perturbation approach, the effective pseudo-spin 
models of the degenerate ground states are obtained. Prom these models, we get the 
energy splitting of the ground states, of which the results are consistent with those 
from exact diagonalization method. 
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I. INTRODUCTION 



In quantum mechanics, quantum tunneling effect is a process by which quantum particles 
penetrate barriers, which are forbidden in classical processes^. It is Gamov who pointed 
out that a single a particle can tunnel through a barrier which introduced "macroscopic 
quantum tunneling" (MQT) into physics firstly. Macroscopic quantum tunneling effects 
have been widely applied to different research fields, such as quantum oscillations between 
two degenerate wells of NH 3 , quantum coherence in one dimension charge density waves, 
macroscopic quantum tunneling effect in ferromagnetic single domain magnets and quantum 
tunneling phenomena in biased Josephson junctions. In general, to find MQT in a system, 
there must exist two or more separated "classical" states with macroscopically distinct. As 
shown in Fig.l, a quantum particle may take a short cut from one well to the other without 
climbing the barrier. 

In this paper we will study a new class of MQT - the MQT in Z 2 topological order. At the 
beginning we give a brief introduction to Z 2 topological order. Topological order is a new 
type of quantum orders beyond Landau's symmetry breaking paradigm^^^, of which there 
are four universal properties : 1) All excitations have mass gap; 2) The quantum degeneracy 
of the ground states depends on the genius of the manifold of the background; 3) There 



FIG. 1: The scheme of a typical macroscopic quantum tunneling process. 

are (closed) string net condensations; 4) Quasi-particles have exotic statistics. All these 
properties are robust against perturbations. Z 2 topological order is the simplest topological 
ordered state with three types of quasiparticles: Z 2 charge, Z 2 vortex, and fermions^. Z 2 
charge and Z 2 vortex are all bosons with mutual ir statistics between them. The fermions can 
be regarded as bound states of a Z 2 charge and a Z 2 vortex. In last ten years, several exactly 
solvable spin models with Z 2 topological orders were found, such as the Kitaev toric-code 
modet^, the Wen-plaquette model^ and the Kitaev model on honeycomb lattice^. 

A decade ago, Kitaev pointed out that the degenerate ground states of a Z 2 topological 
order make up a protected code subspace (the so-called toric-code) free from error-^. In 
Ref.-, topological qubit based on the degenerate ground states of a Z 2 topological order has 
been designed. Then one can manipulate the degenerate ground states by braiding anyons, 
which has becomes a hot issue recently^^^^^^^^^^^ 1 ^ 1 ^. Recently, an alternative 
approach to design TQC is proposed by manipulating the protected code subspace^'ii. The 
key point to manipulate the degenerate ground states is to tune their MQT effect. Thus it 
becomes an interesting issue to study the MQT in Z 2 topological order. 

In this paper, by using a high-order degenerate perturbative approach, we study the 
MQT of the degenerate ground states of Z 2 topological order, taking the Wen-plaquette 
Model as an example. The remainder of the paper is organized as follows. In Sec. II, the 
degenerate ground states of the Wen-plaquette model is classified by the topological closed 
string operators. In Sec. Ill, the dynamics of quasi-particles are studied. In Sec. IV, the 
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MQT of the degenerate ground states of the Z 2 topological order are formalized on a torus 
of different lattices. The numerical results are given to compare with the theoretical results. 
Finally, the conclusions are given in Sec. V. 

II. THE DEGENERATE GROUND STATES AND ITS REPRESENTATION OF 

STRING OPERATORS 

In this section, we study the degenerate ground states of the Wen-plaquette model. The 
Hamiltonian of the Wen-plaquette model is given by 



and g > 0. erf, of are Pauli matrices on site i. The ground states of the Wen-plaquette model 
denoted by Fj, = +1 at each site are known to be an example of Z 2 topological state. The 
ground state energy becomes Eq = —gN where N is the total lattice number-^^. 

In the topological order of the Wen-plaquette model, there exist three types of open string 
operators W C (C), W V (C), Wf(C) corresponding to three types of quasi-particles: Z 2 charge, 
Z 2 vortex, and fermion, respectively^. Here C is a (closed or open) loop. To create a Z 2 
vortex (charge) excitation, one may draw a string state that connects nearest neighboring 
even (odd) plaquettes W V (C) (or W C {C)). Such a string state is created by the following 
string operator Ylc a T where the product Ylc is over all the sites on the string along a loop 
C connecting even-plaquettes (or odd-plaquettes), cjj = y if i is even and aj = x if i is odd. 
For a fermionic excitation, the string operator has a form as Wf(C) = Yln* 7 ^ w hh a string 
C connecting the mid-points of the neighboring links, and i n are sites on the string. l m = z 
if the string does not turn at site i m . l m = x or y if the string makes a turn at site i m . l m — y 
if the turn forms a upper-right or lower-left corner. l m — x if the turn forms a lower-right or 
upper-left corner. It is obvious that the fermionic string can be regarded as a bound state 
of strings of the Z 2 charges and the Z 2 vortices, that is Wf(C) = W C {C)W V {C). If C are 
closed loops, we get condensed closed-string operators of the ground states |^o) as 
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FIG. 2: The topological closed string operators on a torus. The dots denote the crosses of different 
types of strings. 

One can see the detailed definition of the string operators in Ref.-. 

To classify the degeneracy of the ground states, we define three types of topological closed- 
string operators W C (C), W V (C) and Wf(C) = W C (C)W V (C), with C denoting topological closed 
loops. The word 'topological' means that the 'big' loops C surround the torus globally (See 
Fig. 2). One can easily check the commutation relations between the topological closed string 
operators and the Hamiltonian 

[H, W C (C)} = [H, W V (C)} = [H, W f {C)\ = 0. (4) 

For the ground states on a torus of an even-by-even (e*e) lattice, we can define four types 
of elementary topological closed string operators, W v (Cx), W v (Cy), Wf(Cx) and Wf(Cy)- 
Here Cx denotes a closed-loop around the torus along e x -direction and Cy denotes a closed 
loop around the torus along ey-direction. Due to the commutation (or anti-commutation) 
relations between them 

[W v (Cx),W f (C x )} = 0, [W v (Cy),W f (Cy)}=0, (5) 
[W v (Cx),W v (C Y )\ = 0, \W f (Cx),Wf(C Y )] = 0, 
{W v (C x ),W f (Cy)} = 0, {W v (C Y ),W f (Cx)} = 0, 

we may identify W v (Cx), W v (Cy), Wf(Cx) and Wf{C Y ) by pseudo-spin operators rf , rf r|, 
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Pesudo-spin operators 


Cx Cy C-xy 


Z 2 -vortex 
Z 2 -charge 
Fermion 


rf <g) 1 1 <g rf rf ® rf 
rf (g) r| rf ® rf rf <g> r| 
1 (g rf rf ®1 rf <g> rf 



TABLE I: Pseudo-spin representation of the topological closed string operators on an even-by-even 
lattice. 

and rf as 

W V (C X ) = rf ® 1, W V (C Y ) 1 ® rf , (6) 
W>(Cx) ->• 1 ® rf , W>(Cy) -> rf ® 1. 

Thus other five topological closed string operators VF c (Cx), Wc(Cy), W c (Cxy), W v (Cxy) and 
Wf(C X Y) are denoted by rf ® rf , rf ® rf , rf <g> r|, rf ® rf and rf <g> rf , respectively, 

W^ c (Cx) - rf <g> rf , W c (Cy) - rf ® rf , (7) 
W^ c (Cxy) -> rf ® r|, W w (C X y) -> rf ® rf , 
W)(C*y) -> rf ®rf. 

Here Cxy is a closed loop around the torus along diagonal directions. In the table. (I), the 
pseudo-spin representation of the topological closed string operators are illustrated. 

Then as the eigenstates of rf (I = 1,2), the four degenerate ground states are denoted 
by I mi, 7712) =| mi)® | m 2 ). For m/ = 0, we have 

rf I mi) =| m { ), (8) 

and for mi = 1 we have 

rf I m z ) = - I mj). (9) 

Physically, the topological degeneracy arises from presence or the absence of ir flux of fermion 
through the hole. The values of m/ reflect the presence (mi = 1) or the absence (m; = 0) of 
the 7r flux in the hole. 

For the degenerate ground states on an even-by-odd (e * o) lattice, the situation changes. 
Because a Z 2 vortex or Z 2 charge has to move even steps to go back to the same plaquette 
around a torus, we cannot well define a topological closed string operator of Z 2 vortex or Z 2 
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Pesudo-spin operators 


Cx Cy Cxy 


Z2-vortex 
Z 2 -charge 
Fermion 


rf — rf 
rf — rf 
1 rf rf 



TABLE II: Pseudo-spin representation of the topological closed string operators on an even-by-odd 
lattice. 

charge along e^-direction, of which the loop consists of odd number plaquettes. So we can 
only define topological closed string operator of Z 2 vortex and Z 2 charge along e x -direction 
W(Cx) (W(Cx) = W v (Cx) = WdCx)) and the corresponding fermionic string operator along 
e^-direction Wf(C Y ). Due to the anti-commutation relations between W(C X ) and Wf(C Y ), 

{W(C X ), W f (C Y )} = 0, (10) 

we may represent W(C X ) and Wf{C Y ) by pseudo-spin operators rf and rf , respectively 

W{C x )^Tf, W)(Cy)-7?. (11) 

Therefore there are two degenerate ground states | mi) that are the eigenstates of rf . In 
the table. (II), the pseudo-spin representation of the topological closed string operators on 
e* o lattice are shown. 

Similarly, for the degenerate ground states on an odd-by-even (o * e) lattice there are also 
two types of closed string operators, W(Cy) (W(C Y ) = W V (C Y ) = W C (C Y )) and Wf(C x ), 
which can be described by pseudo-spin operators rf and rf. Therefore, the two degenerate 
ground states on an o * e lattice are denoted by | m 2 ) which are the eigenstates of r|. 

For the degenerate ground states on an odd-by-odd (o * o) lattice, since the total lattice 
number is odd, we cannot well define Z 2 vortex or Z 2 charge globally any more. Instead, 
we can only define a mixed topological closed-string operator, W(Cxy) = Ylc a T wnere t ne 
product Ylc i s over a ^ ^ ne sites on the string along a diagonal loop C connecting plaquettes. 
The index Oj = x or y is determined by the position of the plaquettes. Because W(Cxy) 
anti-commutes with Wf(Cx) and Wf(Cy), 

{W(C XY ), W f (C x )} = 0, {W(C XY ), W f (C Y )} = 0, (12) 

we may represent W{Cxy) and Wf(Cx) (or Wf(C Y )) by pseudo-spin operators r x and r z , 



7 



Pesudo-spin operators 


Cx Cy Cxy 


Z2-vortex (/^-charge) 
Fermion 


- - T* 



TABLE III: Pseudo-spin representation of the topological closed string operators on an odd-by-odd 
lattice. 

respectively 

W(C XY ) -H- T X , 
Wf(C X ) = WjiCy) - T Z . 

It is noted that 

W f (C XY ) = W f {C x )Wf{C Y ) = 1. (13) 

Thus the two degenerate ground states on an o * o lattice | m) are the eigenstates of r z . In 
the table. (Ill), the pseudo-spin representation of the topological closed string operators on 
o*o lattice are shown. 

As a result, the degeneracy Q of the ground states of the Wen-plaquette model on lattices 
with periodic boundary condition (on a torus) is dependent on the lattice numbers : Q = 4 
on e * e lattice, Q = 2 on other cases (e * o, o * e and o* o lattices)^^ 1 ^ 1 ^ 1 ^. 

III. PROPERTIES OF QUASI-PARTICLES OF THE WEN-PLAQUETTE 

MODEL 

In this section we study the properties of the quasi-particles of the Wen-plaquette model. 
In this model, Z 2 vortex is defined as = — 1 at even sub-plaquette and Z 2 charge is = — 1 
at odd sub-plaquette. The energy gap of Z 2 charge and Z 2 vortex is 2g. The fermions that 
are the bound states of a Z 2 charge and a Z 2 vortex on two neighbor plaquettes have an 
energy gap of 4g. All quasi-particles in such an exactly solvable model have flat bands. The 
energy spectrums are E v = E c = 2g for Z 2 vortex and Z 2 charge, Ef = 4g for fermions, 
respectively. In other words, the quasi-particles cannot move at all. In particular, there exist 
two types of fermions : the fermions on the vertical links and the fermions on the parallel 
links. 



8 




FIG. 3: The hoppings of Z2 vortex, Z2 charge and fermions. The shadow plaquettes, the striped 
plaquettes and the dots on the links represent Z2 vortices, Z<i charges and fermions, respectively. 

Under the perturbation 

H I = h x Y,^ + h z J2 a *> ( 14 ) 

i i 

the quasi-particles (Z2 vortex, Z2 charge and fermion) begin to hop^^^^^^^^. The 
term h x of drives the Z2 vortex, Z2 charge and fermion hopping along diagonal direction 

i 

&x — &y (See Fig. 3(a)). For example, for a Z 2 vortex living at i plaquette Fj = —1, when of 
acts on i + e x site, it hops to % + e x — e y plaquette denoted by F i+Sx x = — 1, 

F l= -1 - Fi = +1, F i+Sx _£ y = +1 F i+gx _ iy = -1. (15) 

A pair of Z2 vortices at i and i + e x — e y plaquettes can be created or annihilated by the 
operation of af, 

F l= +1 -> F t = -1, F i+ e x -e y = +1 -> F i+8x _, y = -1. (16) 
The term h z ^2 o\ drive fermion hopping along e x and e y directions without affecting Z 2 

i 

vortex and Z2 charge : the fermions on the vertical links move along vertical directions and 
the fermions on the parallel links move along parallel directions. With the help of the term 
h x ^2 a i, the two types of fermions are mixed and the fermions may turn round from vertical 

i 

links to parallel links (See Fig. 3(b)). 

A fact is that the topological closed string operators can be considered as quantum tunnel- 
ing processes of virtual quasi-particle moving along the same loops. Let us take the quantum 
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tunneling process of Z 2 vortex as an example : at first a pair of Z 2 vortices are created. One 
Z 2 vortex propagates around the torus driven by operators and annihilates with the other 
Z 2 vortex. Then a string of af is left on the tunneling path, which is just the topological 
closed string operator W V (C). Such a process effectively adds a unit of a 7r-fiux to one hole 
of the torus and changes mi by 1 . 

IV. MACROSCOPIC QUANTUM TUNNELING EFFECTS OF THE 
DEGENERATE GROUND STATES 

It is known that the degenerate ground states of Z 2 topological orders have the same 
energy in the thermodynamic limit. The different ground states can not mix into each other 
through any local fluctuations. However, in a finite system, the degeneracy of the ground 
states can be (partially) removed due to quantum tunneling processes, of which virtual 
quasi-particles move around the torus^i&iQJi In general cases, one will get large energy 
gaps for all quasi-particles and very tiny energy splitting of the degenerate ground states 
AE. Based on such condition, we may ignore high energy excited states and consider only 
the degenerate ground states. Thus in the following parts we only focus on the ground states 
that are a four-level (or two-level) system. 

A. The high-order degenerate perturbation theory 

To solve quantum tunneling problems, people have developed many approaches includ- 
ing the well known WKB (Wentzel, Kramers and Brillouin) method and the instanton 
approach lately. However, based on semi-classical approximation both above approaches 
are not available to the MQT of Z 2 topological order. Instead, in this part, we develop a 
high-order degenerate perturbative approach to calculate the MQT. 

The Hamiltonian of the Wen-plaquette model under the external field has a form as 

H = H + Hj (17) 
in which H = —g ^\ e>"faf + ~ af +i + ~ af + ~ is the unperturbation term, and Hj = h x ^2 a f + 

i 

h z °f is the small perturbation one. For simplicity, we consider the quantum tunneling 

i 
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process between two degenerate ground states | m) and | n), 

| m) <^>| n). (18) 
According to the Gell-Mann-Low theory, we define a transformation operator Uj(0, — oo) 



as 

r-0 



C7j(0,-oo) = Texp(-i I Hi(t')dt') (19) 

J — oo 

where 

H' T {t) = e^/fie - * 2 * *. (20) 

Here T denotes a time order and h = 1. Then the transformation operator C//(0, — oo) in 
Eq. (1191) can be written as 

oo 

f//(0,-oo) |m) = ^C/{ j) (0,-oo) |m) , (21) 

3=0 

where 



f>i 0) (0,-oo) |m) = |m), 

f>i 1} (0,-oo) |m> — —i Hj{t)dt\m) 

J — oo 



if/ |m) 



-So — #o 



/•O 

t/{ 2) (0,-oo) |m> = —i H'^uPiO.-^dtlm) 

J — oo 



oo 

Hi l -^-H! \m) 



Eq — Hq Eq — Hq 

f/j j¥0) (0,-oo) \m) = ( l -^-H I y\m). (22) 

Eq — Hq 

The element of the transformation matrix from the state \m) to |n) becomes 

(n\ Ui(0, — oo) \m) 

and the corresponding energy is obtained as 

E = (n\ HU^O, -oo) \m) = E + 5E (23) 

where Eq is the eigenvalue of the Hamiltonian Hq of \m) . 

For the tunneling process from | m) to | n), a quasi-particle will move around the torus 
that leads to topological closed string operator behind. So in the sum of j, the dominated 
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term is labeled by j — L — 1. L is the length of the loop of a topological string operator 
W v (Ca) where v = v, c or / and A = X, Y or XY. Then considering the tunneling process 
corresponding to W V (C\), we obtain the perturbative energy as 

5E = (n | F/£7 7 (0, -oo) | m) (24) 

oo 

= (n | ifjj^ C/j^-oo) \m) 

3=0 

= (n | if/L^^-oo) |m) 

Now it is noted that the operator if f L r j L ~ 1 ' ) (0, — oo) is proportion to a topological string 
operator W„(Ca). 

Considering all tunneling processes, we may denote the ground state energies as a four- 
by-four matrix (for the four degenerate ground states on e * e lattice) or two-by-two matrix 
(for the two degenerate ground states on e * o, o * e and o* o lattices), 

5E = J2( n I i?/^/ L-1) (0,-oo) \m) . (25) 

m,n 

Finally we can diagonalize the four-by-four or two-by-two matrices and obtain the energy 
splitting. 

B. Macroscopic quantum tunneling effect of the degenerate ground states on o * o 

lattice 

Firstly, we study the MQT of the two degenerate ground states on an L x x L y (L x and 
L y are odd numbers and L x > L y ) lattice. For simplicity, we use |f) and ||) to describe 
the two degenerate ground states | m — 0) and | m = 1), respectively, of which the two 
ground states can be mapped onto quantum states of pseudo-spin f. Under the perturbation, 
Hj = h x ^2 af + h z Y2 of , two types of quantum tunneling processes dominate - the one that 

i i 

Z2 vortex (or Z% charge) propagates around the torus along diagonal direction and the other 
that fermion propagates around the torus along e y -direction. 

For the first tunneling process, a virtual Z2 vortex (or Z2 charge) will run around the 
torus as long as a path with length L that is equal to Lx ^ y . Here £ is the maximum common 
divisor for L x and L y . For example, on a 3 x 3 lattice, we get L = -^p = 3; on a 3 x 5 
lattice, we get L = ^ = 15. 
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FIG. 4: Generation and Hopping of Z 2 vortex. The shadow plaquettes represent Z 2 vortices. 

From Eq. (I25j) . one may obtain the energy splitting AE of the two ground states as 

8E = UP = (T| M-^rH;)^ ||) . (26) 
E — H 

Due to the translation invariance, to calculate ( 1 - Hi) M) = ( J 1 * f, XVf) ||) , we can 
choose site i as the starting point of the tunneling process and get 

Hj)\i) - L x L y (— *of) ||> (27) 



Eq — Hq Eq — Hq 

= L x L y (— 

Eq — Hq 

where is the excited state of two Z 2 vortices (or Z 2 charges) at plaquettes i — e y and 
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i — e x with an energy E + Ag (See Fig. 4(a)). From H = (E + Ag) , we have 

-H I )\i) = L a L v {-^-)\^ i ). 



In next step, one Z 2 vortex (or Z 2 charge) moves one step, we get 



£°f) a u> 



h x h x 

l x lj — )( — ^y^tf?) 

y -4g J{ E -H l+ex eyJl 1 



L x Ly{—)(—)\%). (28) 



^ K ^4<7' 

where is the excited state of two Z 2 vortices (or Z 2 charges) at plaquettes % + e x — 2e y 
and i — e x . See Fig. 4(b). Then step by step, one Z 2 vortex (or Z 2 charge) moves around the 
torus. When the Z 2 vortex (or Z 2 charge) goes back to its starting point and annihilates 
with the other, the original quantum state ||) changes into |f). Finally we get the energy 
splitting 

A = 25E = 2U\ L) = 2(t| Hi{ 1 -^H I ) L °- 1 1|) 

E — H 

= 2 >< L * L ' (i£)£-i = 8L ^^(^) L °- ( 29 ) 

It is noted that Lq — 1 is an even number. 

Because the quantum tunneling process of Z 2 vortex (or Z 2 charge) plays a role of r x on 

—':)" 

we obtain the effective pseudo-spin Hamiltonian due to the contribution of Z 2 vortex (or Z 2 
charges) as 

^e ff = y(|T)ai + ll)(Tl) = ^ (31) 
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FIG. 5: Tunneling path of virtual fermion along e^-direction on an 5 x 5 lattice (The dots on the 
links denote the fermions). 

where J x = A/2±^±. 

For the second tunneling process, a virtual fermion will move around the torus along 
direction e y with length L y (It is noted that due to L x > L y , the length of tunneling path 
along e x direction is longer). See Fig.5. Such a tunneling process changes the quantum 

/'it A ( IT) \ ( \v\ 

states turn into = t z \ \ . The extra sign of the state M) comes 

from the presence of ir flux of fermionic quasi-particles through the holes of the torus. From 
Eqf551 we can get the energy shift of the state ||) as 

5E = f;<I| M-^-HiY II) = ^tS£ (32) 



j=0 



E — H { 



o 



with an even number L y — 1. Through the same approach, we get the energy shift AE of ||) 



is equal to —L x L y /, ^£ v -i ■ Then an energy difference e of the two ground states is obtained 



as 



e = 25E = lQL x L y g( f f) L y. (33) 

$9 



Finally the two-level quantum system of the two degenerate ground states on an o * o 
lattice can be described by a simple effective pseudo-spin Hamiltonian 

?u = y(ir>ai + u><ti)+|(iTXTi-u>ai) (34) 

= J x r x + J z r z 
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where J x = A/2 and J z = s/2. By diagonalizing the effective Hamiltonian matrix, we can 
get the eigenvalues of the two ground states 



2 ' v 2 

The total energy splitting becomes 



E ± = ± x /£y + ( £ -y. (35) 



/(y) 2 + (|) 2 - (36) 



For the Wen-plaquette model under external field along x-direction, the total energy splitting 
AE is reduced into A = 8L x L y g(j^) L ° . On the other hand, for the Wen-plaquette model 
under external field along ^-direction, the total energy splitting AE is s = 16L x L y g(^) Ly . 

C. Macroscopic quantum tunneling effect of the degenerate ground states on e * o 

lattice 

Secondly, we study the MQT of the two degenerate ground states on an L x x L y (L x is 
an even number and L y is an odd number) lattice^. Now we map the two-fold degenerate 
ground states | mi = 0) and | m\ = 1) onto quantum states of the pseudo-spin fi as |f)i and 
| respectively. Under the perturbation, Hj = h x &f + h z a ti there are two types of 

i i 

quantum tunneling processes - virtual Z 2 -vortex (or Z 2 charge) propagating along e x — e y 
directions around the torus and virtual fermion propagating along e y direction around the 
torus. 

For the virtual Z 2 -vortex (or Z 2 charge) propagating along e x — e y directions around the 
torus, the energy splitting A can be obtained by the high-order degenerate-state perturbation 
theory as 

A = 2<T|i£,(— ^^/) L °- 1 |l)i (37) 
Eo — H 

Because the quantum tunneling process of Z 2 vortex (or Z 2 charge) plays a role of rf on the 
quantum states IT>1 as ^ = * IT>1 , we obtain the effective pseudo-spin 

V IT)l / V l)l J 

Hamiltonian due to the contribution of Z 2 vortex (or Z 2 charge) as 

^e ff = y(|T)i(||i + ||)i(T|i) = J.rf (38) 
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FIG. 6: The energy splitting between the two degenerate ground states of the Wen-plaquette model 
in an external field along x-direction {g = 1). Here N * M denotes a N x M lattice. 

where J x = A/2. 

For the tunneling process of fermion propagating around the torus along direction e y , we 
obtain the energy difference e of the two ground states as 

e = 2AE = lQL x L y g(^) L y . (39) 

The length of the tunneling path is L y which is an odd number. Such tunneling process 
plays a role of rf . 

Finally the two-level quantum system of the two degenerate ground states on an e * o 
lattice can be described by 

Tiff = J x t? + Jz-rf (40) 
where J x = A/2 and J z = e/2. The total energy splitting now becomes 

AE = E + - E. = 2^(y) 2 + (|) 2 . (41) 

In Fig. 6 and Fig. 7, we plot the numerical results from the exact diagonalization technique of 
the Wen-plaquette model on different o*o and e*o lattices. Table. (IV) shows the tunneling 
lengths L from the numerical results (the numbers in the brackets are the theoretical 
predictions), which indicate that our theoretical results are consistent with the numerical 
results from exact diagonalization approach. 
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FIG. 7: The energy splitting between the two degenerate ground states of the Wen-plaquette model 
in an external field along ^-direction (g = 1). Here N * M denotes a N X M lattice. 



L 


3*3 2*5 3*4 3*5 


h x 


2.98312 (3) 9.84653 (10) 12.10754 (12) 15.01707 (15) 


h z 


3.06994 (3) 4.83737 (5) 3.03164 (3) 3.01557 (3) 



TABLE IV: The tunneling lengths Lq from the numerical results (the numbers in the brackets 
are the theoretical predictions). h x means the external field along x-direction and h z means the 
external field along z-direction. Here N * M denotes a N x M lattice. 

D. Macroscopic quantum tunneling effect of the degenerate ground states on e * e 

lattice 

Thirdly, we study the MQT of the four degenerate ground states on an L x x L y {L x 
and L y are even numbers with L x > L y ) lattice. We denote the four degenerate ground 
states | mi,m 2 ) = 1 , 0) , |1,0), |0, 1), |1, 1) by the quantum states of pseudo-spin fi and f 2 . 
Under the perturbation, Hi = h x of + h z ^2 o\, there are five types of quantum tunneling 

i i 

processes - virtual Z^- vortex propagating along e x — e y direction around the torus, Z 2 charge 
propagating along e x — e y direction around the torus, and virtual fermion propagating along 
e x , e y , e x — e y direction around the torus, respectively. We will calculate the ground state 
energy splitting from the degenerate perturbation approach one by one. 
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In the first step we study the quantum tunneling process of Z 2 -vortex propagating along 
e x — e y direction around the torus. After such tunneling process, the quantum states 
f\0,0)\ 
I 1,0) 
1 0,1) 

VI 1 - 1 )/ 



turn into 



f\o,o)\ 
I i,o) 

I 0,1) 

VIM)/ 



/lM)\ 

I o,i) 

I i,o) 
VI o,o)/ 



f\o,o)\ 
1 1,0) 

I 0,1) 

VI 1 - 1 )/ 



(42) 



Thus we may use the pseudo-spin operator t*®t% to denote the tunneling process. The 
effective pseudo-spin Hamiltonian due to the contribution of Z 2 vortex is obtained as 



fteS = JxxTi®T% 



(43) 



where J xx = L x L y ,_ /l -i ■ Similar to the results in above section, the length Lq of the 



x y (-4 fl ) 



tunneling path is equal to Lx ^ v where £ is the maximum common divisor for L x and L y . 

In the second step we study the quantum tunneling process of Z 2 -charge propagating 
along e x — e y directions around the torus. We may use the pseudo-spin operator rf ®t% to 
denote this tunneling process, of which the effective pseudo-spin Hamiltonian is obtained as 



ft. 



eff 



(44) 



where J„ 



L x Li 



' fe ?L n °-i and Lq 



'yy ~ ^x^y (_ 4g )L -i -^u — ^ 
In the third step we study the quantum tunneling process of fermion propagating along 
e x and e y directions around the torus, of which the pseudo-spin operators correspond to 
1 ® r| and rf <8> 1, respectively. Then the effective pseudo-spin Hamiltonian due to the 
contribution of the two quantum tunneling processes is obtained as 



ft. 



off 



h* (rf <g>l) + h z 2 (l®r. 



(45) 



where h\ = -8L x L y g(^) L * and ~h\ = -8L x L y g(^) L v . 

In the last step we study the quantum tunneling process of fermion propagating along 
e x — e y direction around the torus, of which the pseudo-spin operator corresponds to rf®r|, 
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FIG. 8: Tunneling path of virtual fermions on an 6 x 6 lattice with 4 corners (The dots on the 
links denotes fermions). 

respectively. Now there are a lot of tunneling pathes with same length 2L . Different tun- 
neling pathes can be labeled by the positions of corners, at which the fermions make a turn 
round from vertical links to parallel links (or parallel links to vertical links). For a path with 
2k corners (k is an positive integer number), the topological closed string operator can be 
written as 

W f (C XY ) = ^^+1-^-1^^+1-^0-2^0-1^0 ( 46 ) 

with the site i = (i x , i y ) and a neighboring site % + 1. See Fig. 8. Along the closed loops, each 
operator <rj corresponds to a corner. Therefore, the number of pathes with 2k corners that 
is equal to the power of h x in AE is obtained as 

Cl °- 1 ~ k\(L -k-iy: {A7) 

It is noted that for any path, there are at least two corners. Then after considering the 
tunneling processes of all possible pathes, the matrix element of rf ®r| is obtained as 

J — e — L L C 1 ^ 2 ^ 2L0 " 2 (48) 
2 (2h*r(h z ) 2L °- 4 

+L * L V U L -1 (_g^2L -l + •- 

(2h*) 2k (h z ) 2L "- 2k 
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(2h x ) 2L ° 

L * L f(_ % )2Lo-l 



— L X Ly 



\{2h x ) 2 + (h 



z\2]L 



{h z ) 



z\2L 



\2L -1 



Finally we derive an effective pseudo-spin Hamiltonian of the four ground states as 



ft cS ~ J xx (rf ®7f ) + J ro (rf®^) + J Z2 (rf<g)T|) + ^ (rf 01) + h\ (l®r 2 2 ) (49) 
which is equal to 



ft 



( J zz + h\ + h z 2 

—J, 



eff 



/if - 







Jxx J% 



tltl 



yy 





Jxx Jyy 

Jzz hj\-\-h?2 

J z 



Jxx J yy 








h\ -h z 2 J 



(50) 



The coefficients of 7^ e ff are given by 



Jxx Jyy -Lx-Ly 



(h 



x\Lo 



Jzz J J xJ J y 



_4#)£o-i' 
[(2h x ) 2 + (^) 2 ] Lo - (h z ) 2Lo 



(51) 



hi = -8L x L y g(^-) L % h* = -8L x L y g(^) L *. 
By diagonalizing the effective Hamiltonian, we get the energies of the ground states as 



E 2 
E 3 
E, 



Jzz - V (M - K) 2 + AJl 



(52) 



Jzz + \f (hi - hi) 2 + AJ 2 X 

Jzz +K + K, 

Jzz h-y h 2 - 



Because the parameter J zz is always much smaller than others as \J ZZ \ <C \J X x\ , \J% 



yy\ ' 



we may simplify ft cS as^ 



ft eS ~ J xx {t x ®t x ) + J TO (rf ®r|) + ^ (r 2 ®!) + (l®r 2 2 



(53) 



and obtain the energies as 



Ex 



(k - h) 2 + u 2 x 



(54) 
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E 2 ~ ^(hf-h* 2 )* + U* x , 
E 3 ~hl + ~h z 2 , 
E A ~ 

Then when the external field increases (h x 7^ and h z ^ 0), the single energy level of the 
initial four degenerate ground states split into four energy levels. 

If we apply the external field along ^-direction, the four energy levels are 

Ei ~ -hi + h z 2 , Ei~h\- h z 2 , (55) 
E 3 ~ hi + h z , E A ~ -hi - h z 2 , 

where hi = —8L x L y g{^-) Lx and h 2 = — 8L x L y g(^) Ly . In the anisotropic limit, L x L y , 



we have 



hi 



h\ 



. In this case, the initial four degenerate ground states split into two 
groups, Ei ~ h 2 , E 2 — —h 2 , E 3 ~ h 2 and E4 ~ —h 2 . In each group, there are two energy 
levels, of which the energy splitting E\ — E 3 = 2h\ is very tiny. In contrast, the energy 
"gap" between the two groups Ei — E 2 = —2h 2 is larger. One can see the energy levels of 
the Wen-plaquette model in external field along z-direction on 2 x 6 lattice (g = 1) in Fig.9. 
In the isotropic case, L x = L y , we have h\ = h z 2 . In this case, the initial four degenerate 
ground states split into Ei = E 2 = 0, E 3 = 2h\ and E4 = —2h\. One can see the energy 
levels of the Wen-plaquette model in external field along ^-direction on 4 x 4 lattice (g — 1) 
in Fig. 10. 

On the other hand, if we apply the external field along x-direction, the four energy levels 
become 

Ei ~ -2J XX , E 2 ~ 2J XX , (56) 
E 3 = E 4 = J zz ~ 0, 

where J xx = L x L y rj^p^i and J zz = — 8L x L y g(^) 2L ° . Now the initial four degenerate 
ground states split into three energy levels. One can see the energy levels of the Wen- 
plaquette model in an external field along x-direction on 4 x 4 lattice (g — 1) in Fig. 11. 
In addition, one may consider the MQT under a more general perturbation 

& = of + /*» £ o* + h z J2 ( 5 7) 
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FIG. 10: The ground state energies of the Wen-plaquette model in an external field along z-direction 
on 4 x 4 lattice (g = 1). 

For an external field of h x ^ 0, h y ^ and h z ^ 0, all quasi-particles (Z2 vortex, Z2 charge 
and fermion) can move along e x , e y , e x ± e y directions freely. Therefore, to calculate the 
MQT of the degenerate ground states on an e * e lattice, all nine types of quantum tunneling 
processes should be considered. The corresponding effective pseudo-spin Hamiltonian of the 
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FIG. 11: The ground state energies of the Wen-plaquette model in an external field along x-direction 
on 4 x 4 lattice (g = 1). 

four ground states turns into 



where J xx , J yy} J zz , J ZX} J xz , hf, /if, hf, h\ are determined by the energy splitting of the 
degenerate ground states from the nine tunneling processes^ 1 ^. This issue (the MQT of Eq. 
( 1571) ) will be studied elsewhere. 



In this paper, we study macroscopic quantum tunneling (MQT) effect of topological 
order in the Wen-Plaquette model that is characterized by the quantum tunneling processes 
of different virtual quasi-particles moving around the torus. By focusing on the degenerate 
ground states, we get their effective pseudo-spin models. The coefficients of these effective 
pseudo-spin models are obtained by a high-order degenerate perturbation approach. With 
the help of the effective pseudo-spin models, the energies of the ground states are calculated 
and the results are consistent with those from exact diagnalization numerical technique. 



Jxx (t? ®t|) + J yy (rf ®r|) + J zz (rf®7^) 
+J ZX (?f ®7f ) + J xz (rf®r|) + hi (rf ®1) 




(58) 



V. CONCLUSION 
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In the future, the approach will be applied onto the MQTs of Z 2 topological order in other 
models, such as the Kitaev toric-code mode and the Kitaev model on honeycomb lattice. 
After learning the nature of the MQT of Z2 topological orders in different models, one may 
know how to manipulate the degenerate ground states by controlling the external field and 
then do topological quantum computation within the degenerate ground states 1 ^ 1 -. 
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